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ON  SYMMETRIC  REINFORCEMENTS  OF 
CIRCULAR  CUTOUTS1 2 3 
by 

2 

E.  Levin 


A slab  with  a circular  cutout  is  subjected 
to  stresses  in  the  plane  of  the  slab.  The 
cutout  has  removed  material  which  would 
have  participated  in  carrying  the  load, 
hence  the  slab  with  cutout  will  fail  under 
the  application  of  stresses  v/hich  the  com- 
plete slab  could  have  supported.  In  order 
to  eliminate  at  least  part  of  tills  weaken- 
ing, the  slab  with  cutout  may  be  reinforced 
by  the  addition  of  material  about  the  cut- 
out. Such  a reinforcement  may,  of  course, 
be  designed  in  any  shape.  The  present 
paper  is  concerned  with  extending  the  re- 
sults of  Weiss,  Prager,  and  Hodge  [1]3  for 
a cylindrical  reinforcement  to  a reinforce- 
ment of  arbitrary  symmetric  cross-soction. 
In  particular,  a method  will  be  described 
for  the  determination  of  a lower  bound  on 
the  collapse  load  for  a slab  with  circular 
cutout  and  a general  symmetric  reinforce- 
ment. 


1*  lQ£E2&fl£fciaB« 

Consider  a uniform  plane  rectangular  slab  with  a centered 
circular  cutout  of  radius  a.  The  cutout  is  reinforced  by  addition- 
al material  extending  to  the  radius  a + 6;  the  cross-section  of 
the  reinforcement  is  prescribed  (Fig,  1),  The  slab  is  subjected 


1.  The  results  presented  in  this  papor  were  obtained  in  tho 
course  of  research  conducted  under  Contract  N7onr-358lO 
between  the  Office  of  Naval  Research  and  Brown  University. 

2.  Graduate  Student,  University  of  California,  Los  Angeles. 

3.  Numbers  in  square  brackots  refer  to  tho  Bibliography  at  the 
end  of  the  papor. 
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to  arbitrary  uniform  tensions  along  its  outer  edges  and  the  inner 
edge  is  stress  free. 

The  present  paper  is  concerned  with  the  determination 
of  a lower  bound  on  those  values  of  the  loads  which  may  be  con- 
sidered safe.  In  the  next  section,  the  physical  assumptions  made 
in  discussing  the  problem  will  be  stated,  and  a precise  mathemati- 
cal problem  formulated.  In  particular,  a domain  of  safe  loads 
will  be  defined.  Section  3 is  concerned  with  the  problem  of  equal 
biaxial  tensions.  The  next  three  sections  are  concerned  with  a 
reinforcement  whose  cross-section  is  symmetric  about  both  axes 
through  its  center  and  whose  boundary  in  the  first  quadrant  is 
monotonically  non-increasing. * Some  general  results  are  estab- 
lished for  all  symmetric  reinforcements  by  approximating  the  true 
interaction  curve  by  a parabola.  It  is  shown  that  the  cylindrical 
reinforcement,  which  was  first  solved  by  Weiss,  Prager,  and  Hodge 
[1],  appears  as  a special  case.  An  estimate  of  the  error  intro- 
duced by  the  parabolic  approximation  is  obtained  fpr  a quasi- 
toroidal  reinforcement.  Limitations  and  conclusions  concerning 
the  results  are  prosented  in  section  7. 

2*  Statement  of  the  Problem. 

The  analysis  will  be  based  on  a limit  theorem  of  Drucker, 
Greenberg,  and  Prager  [2]  which  says  that  a given  system  of  loads 
is  safe  if  stresses  can  be  found  which 

1)  are  in  equilibrium  with  the  given  loads,  and 

2)  do  nowhere  exceed  tho  limits  set  by  tho  yield  criterion# 

* Non- symmetric  reinforcements  may  bo  treated  numerically#  The 
method  is  illustrated  with  a bevelled  reinforcement  in  [53* 


The  slab  and  the  reinforcement  are  assumed  to  be  of  a uniform 
perfectly  plastic  material  which  satisfies  Tresca's  yield  criter- 
ion [3]. 

Following  Weiss,  Prager,  and  Hodge  [1],  the  reinforced 
slab  will  be  regarded  as  consisting  of  a uniform  part  under  con- 
ditions of  plane  stress,  and  a hub  which  behaves  as  a curved 
beam  (Fig.  lc).  Let  the  applied  tractions  on  the  oute**  edge  of 
the  slab  be  X^s  and  \ys,  where  s is  the  yield  strength  of  the 
material  in  simple  tension  and  it  is  assumed  that  the  applied 
loads  would  be  safe  if  there  were  no  cutout.  It  follows  that  in 
the  uniform  part  of  the  slab,  the  homogeneous  state  of  stress 
o = s\x,  o = s\v,  t = 0 satisfies  the  two  conditions  of  the 
Drucker,  Greenberg,  Prager  theorem,  and  transmits  total  loads 
of  Intensities  hsXx,  hsXy  per  unit  length  to  the  hub.  The  prob- 
lem, then,  is  to  find  values  of  \x,  \y  which  lead  to  a safe  stress 
distribution  in  the  hub. 

It  has  previously  been  shown  [l,  4]  that  the  safe  states 
of  loading  for  the  uniform  slab  without  cutout  may  be  represented 
by  the  points  on  and  Inside  a hexagon,  and  that  any  such  safe 
state  can  be  obtained  from  the  stress-free  state,  the  state  of 
uniaxial  tension  and  the  state  of  biaxial  tension,  by  the  follow- 
ing operations: 

(a)  If  Tx,  Ty  is  a safe  load,  then  «TX,  -Ty  is  also  a safe 

load. 

(b)  If  Tx,  Ty  and  Tx,  Ty  are  each  safe  loads,  then 

l*Tx  + (1  - p.)Tx,  yTy  + (1  - ^)Ty  is  a safe  load  for 
any  0 < p < 1. 

For  the  roinforcod  cut  slab,  the  domain  of  safe  loads  will  be 
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defined  as  the  largest  hexagon  which  may  be  obtained  from  the 
stress-free,  uniaxial,  and  biaxial  states  by  these  operations. 
Hence  the  problem  of  determining  the  domain  of  safe  loads  is 
reduced  to  the  problems  of  finding  lower  bounds  for  X in  biaxial 
tension  (Xx  = Xy  = X)  and  uniaxial  tension  (Xx  = X,  Xy  = 0). 

The  case  of  biaxial  tension  is  fairly  simple  and  may  be  treated 
for  all  reinforcements  simultaneously.  Uniaxial  tension  is  more 
difficult  and  only  symmetric  reinforcements  will  be  considered. 

For  a curved  beam  the  state  of  stress  in  any  section 
is  specified  by  a bending  moment  M and  an  axial  force  N caused 
by  tensile  stresses,  and  a shear  force  caused  by  shear  stresses. 
The  effect  of  the  shear  force  is  negligible  for  solid  cross  sec- 
tions and  will  be  disregarded  in  the  subsequent  analysis. 

3.  Bjjyaai  T<nsiQQ- 

Consider  a cutout  slab  shown  in  Fig.  1,  and  let  A be 
the  total  crocs-sectlonal  area  of  the  hub.  In  biaxial  tension 
the  hub  is  subjected  to  a uniformly  distributed  radial  force  of 
magnitude  hXg  per  unit  length  of  She  circumference  of  radius 
(a  + t>).  To  satisfy  the  equilibrium  requirement,  the  stress 
resultants  in  the  hub  must  be  given  by 

N = (a  + b)hXs,  M = 0 (3.1) 

at  all  cross  sections.  Such  resultants  are  furnished  by  a uni- 
form tensile  stress  of  magnitude  N/A.  From  tho  yield  criterion 
this  stress  cannot  exceed  the  yield  stress  s,  hence  the  axial 
forco  must  satisfy 


N < As. 


(3.2) 
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Combining  (3.1)  and  (3.2)  leads  to  the  inequality 

\ < = A/[h(a  + & )].  (3.3) 

Hence  for  the  case  of  biaxial  tension  any  value  of  \ satisfying 
Inequality  (3»3)  furnishes  a lower  bound* 


4.  General  Symmetric  Reinforcement* 

The  type  of  reinforcement  considered  here  has  the  follow- 
ing properties: 

1)  the  hub  is  symmetric  about  the  x and  z axes,  where  the 
x axis  lies  in  the  plane  of  the  slab  and  passes  through  the 
center  of  mass  or  the  hub  and  through  the  center  of  the  cutout, 
and  the  z axis  is  normal  to  the  slab  and  also  passes  through  the 
center  of  mass  of  tho  hub,  and 

2)  in  the  first  quadrant  the  boundary  of  the  hub,  z = f(x), 
is  a continuous,  single-valued,  monotonically  non- increasing 
function  of  x (Fig,  2a). 

To  find  limiting  combinations  of  bending  moment  M and 
axial  force  N which  the  beam  can  support,  i.  e. , which  do  not 
violate  the  yield  criterion,  consider  a fully  plastic  stress 
distribution  of  the  type  shown  in  Fig.  2b.  Part  of  the  section 
is  stressed  to  the  yield  limit  in  tension  and  the  rest  in  com- 
pression. The  absolute  values  of  the  bending  moment  and  axial 

force  for  such  a stross  distribution  are 

6/2  r r 

I N|  = | s( 2 ; f (x)dx  - I 2 f (x)dx)|  = |4sl  f(x)dx,  (4,1) 

J £ -5/2  uO 


,V2 

1 

.f(x)  f 

if(x) 

* 

1 Ml  = | s(2  | 

xdxd?  - 2 ' 

i 

xdxdz)|  = 1 4s 

0 6/2 ^ 

0 t 

,6/2 

x f(x)dx|. 
C (4,2) 
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If  N/s  and  M/s  are  regarded  as  Cartesian  co-ordinates 
in  a stress-resultant  space,  then  (4.1)  and  (4.2)  define  an  inter- 
action curve  in  the  NM  plane.  It  can  be  shown  that  this  curve 
will  be  closed  and  will  bound  a convex,  symmetric  region  which 
will  be  defined  as  the  safe  region.  Any  point  within  the  safe 
region  represents  a bending  moment  and  an  axial  force  which  do 
not  violate  the  yield  criterion.  The  boundary  of  the  region  is 
described  parametrically  in  terms  of  C and  represents  the  limit- 
ing admissible  combinations  of  bending  moment  and  axial  force. 

The  exact  safe  region  is  defined  by  Eqs.  (4.1)  and 
(4.2),  which  contain  integrals  of  the  arbitrary  function  f(x)  and 
xf(x).  Therefore,  it  is  not  possible  to  obtain  a direct  limiting 
relation  between  N and  M,  However,  any  region  which  lies  wholly 
within  the  safe  region  will  itself  consist  only  of  safe  stress 
resultants.  It  will  be  shown  that  the  domain  bounded  in  each 
quadrant  by  a certain  parabolic  arc  is  such  a safe  approximation. 

From  symmetry,  it  is  only  necessary  to  consider  the 
first  quadrant.  Here 


Define 


M/s 


N/s  = 4 


b/2 

[ 

= 4 xf(x)dx, 
£ 


f(x)dx. 


Vs 


Nq/s 


= 4 


= 4 


xf (x)dx, 


0 

h(? 


f (x)dx, 


(4.3) 


(4.4) 


(4.  5) 


(4.6) 
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where  M0  and  NQ  represent  the  extreme  values  in  the  positive  M 
and  N directions  respectively.  Thus  the  end  points  of  the  inter- 
action curve  in  the  first  quadrant  are  (0,  M0/s)  and  (No/s,0) 
which  correspond  to  C = 0 and  C = 6/2  respectively.  The  approxi- 
mating curve  will  be  taken  as  the  parabola 

M/M0  = 1 - ('VN0)2  (4.7) 


which  passes  through  the  same  endpoints  with  a horizontal  tangent 
at  N = 0, 

In  order  to  prove  that  this  curve  is  an  admissible 
approximation  it  is  sufficient  to  show  that  for  any  given  value 
of  C between  zero  and  6/2,  the  ordinate  M/M0  of  the  approximate 
curve  never  exceeds  the  ordinate  M/M0  of  the  interaction  curve. 

In  terms  of  C,  the  ordinate  of  the  approximate  curve  is 

r j*  f(x)dx*,'2 

M/M0  = 1 - (N/N0)  = 1 - 100 


, 575 

I J*  f(x)dx 


(4.8) 


while  the  ordinate  of  the  interaction  curve  is 

t*/2 


M/M-  = *5- 


xf (x)dx 


o " y* 


xf (x)dx 


(4.9) 


Therefore  for  all  0 < C < 6/2,  it  must  be  shown  that 


f\6/2 


f xf (x)dx 

e(V  j'4 


r 

JO 


- 1 + 


xf (x)dx 


f 

Jo 


f (x)dx 


-,2 


6/2 

f f(x)dx 
JO  J 


2 0. 


(4. 10) 


It  is  evident  from  the  definition  of  the  approximate 
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curve  that  g(0)  = g(^)  = 0,  and  g * (0)  = 0.  In  the  appendix  it 
is  shown  that  g*(i)  is  equal  to  or  less  than  zero,  and  that  the 
equation  g'(C)  = 0 has  not  more  than  two  roots  in  the  closed 
interval  0 < 6/2.  Since  g(£)  is  continuously  differentiable, 

it  '’ollows  that  in  addition  to  its  extremum  at  C = 0,  it  has 
only  one  relative  extremum  between  zero  and  6/2.  Since  the  slope 
is  non-positive  at  C = b/2,  this  extremum  must  be  a maximum. 
Therefore,  Inequality  (4,10)  is  valid,  since  otherwise  g ( C ) 
would  have  a minimum  in  the  interior. 

Thus  the  approximation  is  safe  and  the  yield  condition 
for  the  hub  may  be  approximated  by 

| M/Mq|  < 1 - (N/Nq)2.  (4.11) 

Under  uniaxial  tension  (Fig.  3)  the  stress  resultants 
M and  N at  any  cross  section  © of  the  hub  may  be  computed  from 
equilibrium  considerations  to  within  a redundant  bending  moment, 
X.  This  statically  indeterminate  quantity  represents  the  bending 
moment  at  the  cross  section  © = 0 and  from  equilibrium  is  seen 
to  be  the  only  non- vanishing  stress  resultant  thero.  The  stress 
resultants  arc  found  to  be 


N = F sin  © = h\s(a+5) sind©, 
fl(a+b)sin© 


(4.12) 


M = X + 


'JO 


h\s  [(a+6/2)sin  0-y]dy=X+  1 h\sa(a+b)sin  0, 

2 (4.13) 


N/N0  = d sin2©,  (4.14) 

M/Mq  = Y + c sin2©, 


(4.1?) 
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where  Y * X/M0,  c(\)  = h\sa(a+6)/2M0  > 0,  d(\)  = h\a(a+6)/N0  > 0* 
The  substitution  of  Eqs,  (4.14),  (4,15)  into  the  yield  condition 
(4,11)  leads  to 

Y + c sin2©  < 1 - d2sinlf©t  (4.16a) 
-Y  - c sin2©  < 1 - d2sin4©.  (4.16b) 


These  inequalities  must  be  satisfied  for  all  0 < © ^ n/2, 

2 

Let  w = sin  ©•  Then  (4, 16)  may  be  written 
' P(w)  = d2w2  «►  cw  ♦ Y - 1 < 0,  (4,17a) 

Q(w)  = dV  - cw  - Y - 1 < 0.  (4,17b) 


Here  P and  Q represent  parabolas  concave  upwards,  hence  a necess- 
ary and  sufficient  condition  that  Inequalities  (4,17)  hold  for 
all  0 < w < 1,  is  that  they  hold  at  both  end  points.  Thus 


- 1 < Y (4,18) 

d2  - c - 1 < Y (4.19) 

Y < 1 (4,20) 

Y < 1 - c - d2.  (4,21) 


The  actual  choice  of  Y is  innaterial,  provided  only 
that  some  Y exists  which  satisfies  all  of  the  above  inequalities, 
A necessary  and  sufficient  condition  that  this  be  the  case  is 
that  each  left  hand  sldo  be  less  than  each  right  hand  side. 

Thus 
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(4.  22) 


- 1 

< 1 - c - d2, 

(4.23) 

d2  - c - 1 

< 1 - c - d2 , 

(4. 24) 

d2  - c - 1 

< !• 

(4.25) 

Since  c and  d are  positive,  the  above  inequalities  are  equivalent 
to 


d + c <2, 


(4.26) 


d < 1. 


(4.27) 


In  view  of  the  definition  of  d,  (4.27)  leads  to 

H2 


\ < 4 


f (x)dx/h(a  + b ) s X. 


(4.28) 


Jo 


where  X^  is  defined  in  Eq.  (3.3). 

With  the  aid  of  the  definitions  of  d,  c,  and  X^,  to- 
gether with  D = 4M0/Nq,  (4.26)  may  be  written 


F(X ) = X2  + (2a/D)X1X  - 2X2  < 0. 


(4.29) 


F(X)  represents  a parabola  concave  upwards,  hence  (4,29)  will  be 
satisfied  for  all  X between  the  two  roots  of  F(X)  = 0.  Since 
the  smaller  root  is  negative  and  X is  positive,  the  restriction 
on  X becomes 

X < \2  = \[  ^2  + (a/D)2  - (a/D)].  (4.30) 


Thus  any  value  of  X satisfying  both  (4,28)  and  (4.30) 
furnishes  a lower  bound  for  uniaxial  loading.  It  follows  from 
(4.28)  and  ( 3*  3 > that  if  X is  safe  for  uniaxial  loading  it  is 
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also  safe  for  biaxial  loading.  Therefore,  the  domain  of  safe 
loads  is  fullv  defined  by 


\ < min(Xi»  X2), 


(4.31) 


where  X^  is  g'.ven  by  (3*3)  and  ^2  ^ (4.30).  In  particular  it 
may  readily  b?  shown  that  (4.31)  is  equivalent  to 


X < X 

x <x2 


if 

if 


a/D  < 1, 
a/D  > 1. 


f\ 


(4.32) 


5.  Cy  1,1  Qlrjc  aj  , Iji 03  Rein £g££&lgQ$ » 

The  results  of  the  preceding  section  are  readily  appli- 
cable to  the  case  of  a reinforcement  consisting  of  a hollow 
cylindrical  ring  about  the  cutout  (Fig.  4).  The  hub  cross  sec- 
tion is  now  a rectangle  and  * = f(x)  is  a constant,  H/2,  where  H 
represents  the  total  height  of  the  reinforcement.  Hence 


X H* 

1 h(a  + 6) 


D = 5, 


H6 


x irr^ — x if  f b 

h(a  +6)  6 ^ 


(5.1) 

(5.2) 

(5.3) 


x ^ <r  viri7S7) 


1 + 2(  6/a)‘ 


-iT 


L 


if  •§  (5.4) 


X = 


For  a "full  strongth"  reinforcement  of  this  type, 
1,  and  the  ring  must  be  so  designed  that 


H/h  = If  (a/6)  ^ 


Il/h  = 


1 + (6/a) 


\fl  + 2(6/a)^  - 1 


iff 


(5.5) 

(5.6) 
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These  results  coincide  with  those  obtained  directly  by  Weiss, 
Prager,  and  Hodge  [1].  The  exact  agreement  is  accounted  for  by 
the  fact  that  for  this  case  only  the  actual  safe  region  in  the 
stress  resultant  plane  is  bounded  by  an  interacting  curve  which 
consists  of  parabolic  arcs.  Hence  the  parabolic  approximation 
used  in  section  4,  coincides  with  the  true  interaction  curve. 

6.  Qua_sl-_t;oroldal  Reinforcement. 

Consider  a reinforcement  of  the  type  shown  in  Fig.  5 
where  the  boundary  of  the  hub  cross  section  is  an  arc  of  a circle, 
f(x)  = V r2  - x2.  A lower  bound  may  readily  be  obtained  from 
Inequalities  (4.28)  and  (4.30).  Further,  an  upper  bound  may  be 
obtained  for  the  error  introduced  by  approximating  the  true 
interaction  curve  by  a parabolic  arc. 

It  can  be  shown  that  the  interaction  curve  in  the  first 
quadrant  is  given  by 

\ u = cos'1(v1)1/3  + (v1)1/3  Vl  - (vx)2/3',  (6.1) 

p 

where  u = N/sr 

and  ^ M/sr3  + (h/2r)^. 

The  appropriate  parabolic  approximation  is 

2 

(v/v0)  = 1 - (u/u0)  (6.2) 

whero 

v = M/sr3  = ^ [ v-^  - (h/2r)3], 
v0  = i Cl  - (h/2r>3), 

uo  = V®*  = 2 [cos_1(h/2r)  + (h/2r)  \fl  - (h/2r)2']. 
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The  curve  (6.1)  may  be  plotted  in  a uv^  space.  However, 
the  parabolic  approximation  as  well  as  the  actual  position  of  the 
u axis  depends  upon  the  ratio  h/2r.  To  obtain  an  estimate  of  the 
error  made  by  using  the  parabolic  aporoximation  (6.2),  a second 
parabola  is  passed  through  the  points  (u  = 0,  v-^  = 1),  (u  = n, 
v^  = 0).  Such  a parabola  is 

~ 1 - (u/x)2.  (6.3) 

It  can  be  shown  by  methods  similar  to  those  employed  in 
section  4 that  the  limiting  parabola  (6.3)  lies  below  the  approx- 
imating parabola  (6.2)  for  all  values  of  h/2r,  so  that,  independ- 
ently of  the  dimensions  of  the  reinforcement,  the  approximating 
parabola  must  lie  between  the  limiting  parabola  and  the  interac- 
tion curve.  As  may  be  seen  from  Fig.  6 or  Table  I the  approxima- 
tion is  a close  one.  The  dotted  curve  represents  a typical 
approximating  parabola  for  the  case  h/2r  = .5* 

- At  least  for  this  type  of  reinforcement  (and  of  course 
for  the  cylindrical  rings  of  section  5)  the  approximation  appoars 
to  be  ve -y  satisfactory.  In  view  of  the  rather  drastic  physical 
approximations  vhich  were  introduced  in  sotting  up  the  problem 
in  section  2,  it  appears  entirely  reasonable  to  adopt  the  para- 
bolic aiproximaMon  to  the  yield  condition. 

7 . Co.ns.Ius ions  and  Limitations. 

The  results  in  section  4 for  a symmetric  reinforcement 
provide  methods  for  obtaining  an  approximate  lower  bound  on  the 
collapse  load.  The  principal  approximation  made  was  that  of 
treating  the  hub  as  a curvod  beam  with  no  shear.  Although  this 
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has  proved  successful  in  the  elastic  range  [6]  certain  limitations 
must  be  observed  when  applying  the  results.  If  the  maximum 
thickness  of  the  hub  is  much  larger  than  that  of  the  slab,  the 
question  of  the  carrying  capacity  arises.  If  frhe  other  hand, 
the  width  becomes  too  large  in  comparison  with  the  inner  radius, 
the  entire  concept  of  a curved  beam  becomes  questionable.  Buck- 
ling in  compression  is  not  considered  and  must  be  treated  sep- 
arately. 

The  above  approximations  tend  to  predict  an  optimistic 
load  capacity.  On  the  other  hand,  the  assumption  of  perfect 
plasticity  which  neglects  strain-hardening  will  make  the  estimate 
conservative.  Also,  the  slab  was  regarded  as  being  in  a homo- 
geneous stato  of  stress  and  no  stress  concentrations  in  the  slab 
wore  considered. 

At  present  there  is  insufficient  experimental  evidence 
available  to  determino  the  extent  to  which  these  various  in- 
fluences tend  to  cancel* 


Appendix. 

The  differentiation  of  Eq.  (4.10)  with  respect  to  C 

loads  to 


- Cf(  C) 

g'(t)  = + 2 


f (x)dx 


0 

and  „5 , f(|)32s2 

e(2)=-vT- 


* I w&/2  2 

xf(x)dx  ( f(x)dx) 

Lj  0 


f(C) 


(Al) 


r 6/2  6 

\ 

(2  i 

j xf(x)dx  - ^ 

f (x)dx 

-Jo  v 

0 

( A2) 


The  sign  of  g'(&)  will  bo  the  same  as  the  sign  of  tho  bracketed 


Bll-17 


15 


expression  since  the  factor  outside  the  brackets  is  always 
positive.  By  combining  the  two  integrals  and  performing  a change 
of  variable,  t = (x  - 6A),  it  is  found  that 

2 1 6/4 

g ’ (6/2)  = j ,P  tJ'f(&A  + t)  - f(5A  - t)|  dt.  (A3) 

Mo  No  Lj0  1 J 

But  throughout  the  range  of  integration  f(&A  + t)  < f(&A  - t) 
since  f is  monotonically  non-increasing.  Hence  the  integrand  is 
non-nositive  throughout  the  range  of  integration  and 


g‘(6/2)  < 0. 


(A*f) 


In  fact  g'(6/2)  = 0 can  hold  if  and  only  if  f(&A  + t)  = f(6/4  -t) 
for  all  0 < t < 6/4,  i.  e.  , f(C)  = c and  the  hub  is  a rectangle. 

For  this  case  the  interaction  curve  is  itself  parabolic  and 


g(C)  = 0. 


Consider  the  number  of  solutions  of  g'(C)  = 0.  Since 


f (0  ? 0,  the  problem  is  to  determine  the  number  of  solutions  of 

5. 


h(C)  = KC  - 


f(x)dx  = 0, 


(A5) 


where  K = [ 


6/2 


f(x)dxf/2[ 


6/2 


xf(x)dx].The  number  of  distinct 


0 


zeros  of  h(C)  may  exceed  by  one  the  number  of  sign  changes  of  its 
derivative 


hf(0  = K - f(C). 


(A6) 


Since  f is  monotonically  non-increasing  and  it  is  constant, 
h’(C)  ^y  change  si  yn  only  once.  Therefore  h(C)  and  hence  g'(C) 
have  at  most  two  distinct  zeros.  3ut  g'(0)  = 0,  thus  there  is 
only  one  extremum  in  the  interval  which  must  represent  a maximum 
since  g * ( 6/2 ) 0 whenever  g(£  ) ^ 0. 
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